General optimality of the Heisenberg limit for quantum metrology 
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Quantum metrology promises improved sensitivity in parameter estimation over classical procedures. However, 
there is an extensive debate over the question how the sensitivity scales with the resources (such as the average 
photon number) and number of queries that are used in estimation procedures. Here, we reconcile the physical 
definition of the relevant resources used in parameter estimation with the information-theoretical scaling in 
terms of the query complexity of a quantum network. This leads to a completely general optimality proof of the 
Heisenberg limit for quantum metrology. We give an example how our proof resolves paradoxes that suggest 
sensitivities beyond the Heisenberg limit, and we show that the Heisenberg limit is an information-theoretic 
interpretation of the Margolus-Levitin bound, rather than Heisenberg's uncertainty relation. 
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Parameter estimation is a fundamental pillar of science and 
technology, and improved measurement techniques for pa- 
rameter estimation have often led to scientific breakthroughs 
and technological advancement. Caves [1] showed that quan- 
tum mechanical systems can in principle produce greater sen- 
sitivity over classical methods, and many quantum parame- 
ter estimation protocols have been proposed since [2]. The 
field of quantum metrology started with the work of Helstrom 
[3, 4], who derived the minimum value for the mean square 
error in a parameter in terms of the density matrix of the 
quantum system and a measurement procedure. This was a 
generalisation of a known result in classical parameter esti- 
mation, called the Cramer-Rao bound. Braunstein and Caves 
[5] showed how this bound can be formulated for the most 
general state preparation and measurement procedures. While 
it is generally a hard problem to show that the Cramer-Rao 
bound can be attained in a given setup, at least it gives an up- 
per limit to the precision of quantum parameter estimation. 

The quantum Cramer-Rao bound is typically formulated 
in terms of the Fisher information, an abstract quantity that 
measures the maximum information about a parameter <p that 
can be extracted from a given measurement procedure. One 
of the central questions in quantum metrology is how the 
Fisher information scales with the physical resources used in 
the measurement procedure. We usually consider two scal- 
ing regimes: First, in the standard quantum limit (SQL) [6] or 
shot-noise limit the Fisher information is constant, and the er- 
ror scales with the inverse square root of the number of times 
T we make a measurement. Second, in the Heisenberg limit 
[7] the error is bounded by the inverse of the physical re- 
sources. Typically, these are expressed in terms of the size N 
of the probe system, e.g., (average) photon number. However, 
it has been clearly demonstrated that this form of the limit 
is not universally valid. For example, Beltran and Luis [8] 
showed that the use of classical optical nonlinearities can lead 
to an error with average photon number scaling N~^l 2 . Boixo 
et al. [9] devised a parameter estimation procedure that sees 
the error scale with N~ k with k € N, and Roy and Braunstein 
[10] construct a procedure that achieves an error that scales 
with 2~ N . The central question is then: What is the real fun- 



damental Heisenberg limit for quantum metrology? We could 
redefine this limit accordingly to scale as 2~ N , but in practice 
this bound will never be tight. 

In this Letter, we give a natural definition of the relevant 
physical resources for quantum metrology based on the gen- 
eral description of a parameter estimation procedure, and we 
prove the fundamental bound on the mean squared error based 
on this resource count. We will show that the resource count 
is proportional to the size of the probe system only if the in- 
teraction between the object and the probe is non-entangling 
over the systems constituting the probe. First, we study the 
query complexity of quantum metrology networks, which will 
lead to a resource count given by the expectation value of 
the generator of translations in the parameter (p. Second, we 
prove that the mean error in <p is bounded by the inverse of 
this resource count. We argue that this is the fundamental 
Heisenberg limit for quantum metrology. Furthermore, we 
show that it is a form of the Margolus-Levitin bound, as op- 
posed to Heisenberg's uncertainty relation. Finally, we illus- 
trate how this general principle can resolve paradoxical situa- 
tions in which the Heisenberg limit seems to be surpassed. 

The most general parameter estimation procedure is shown 
in Fig. la). Consider a probe system prepared in an ini- 
tial quantum state p(0) that is evolved to a state p(<p) by 
U(q>) = e\p(~i(pJif). This is a unitary evolution when we in- 
clude the relevant environment into our description, and it in- 
cludes feed-forward procedures. The Hermitian operator 
is the generator of translations in <p, the parameter we wish 
to estimate. The system is subjected to a generalized mea- 
surement M, described by a Positive Operator Valued Mea- 
sure (POVM) that consists of elements E x , where x denotes 
the measurement outcome. These can be discrete or con- 
tinuous (or a mixture of both). The probability distribution 
that describes the measurement data is given by the Born rule 
p(x\q>) = Ti[E x p((p)}, and the maximum amount of informa- 
tion about <p that can be extracted from this measurement is 
given by the Fisher information 
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FIG. 1: a) General parameter estimation procedure involving state 
preparation P, evolution U (<p) and generalized measurement M with 
outcomes x, which produces a probability distribution p(x\<p). In 
terms of quantum networks, the evolution can be written as a number 
of queries of the parameter <p. b) Example for A' = 4 of the usual 
situation described by -j#glm- where each system performs a single 
query, and the number of queries equals the number of systems (the 
grey box represents Oj((p)); c) for ,-^bfcg the number of queries Q 
does not always equal the number of systems: any two systems can 
jointly perform a single query, and the number of queries then scales 
quadratically with the number of systems; d) for ,j?rb all possible 
subsets of systems perform a single query. The number of queries 
scales exponentially with the number of systems. 



This leads to the quantum Cramer-Rao bound [3,5] 
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where (<S<p) 2 is the mean square error in the parameter <f>, 
and T is the number of times the procedure is repeated. The 
SQL is obtained when the Fisher information is a constant 
with respect to T, and the Heisenberg limit is obtained in a 
single-shot experiment (T = 1) when the Fisher information 
scales quadratically with the resource count. The SQL and the 
Heisenberg limit therefore relate to two fundamentally differ- 
ent quantities, T and F, respectively. We need to reconcile the 
meaning of these two limits if we want to compare them in a 
meaningful way. 

To solve this problem, we can define an unambiguous re- 
source count for parameter estimation by recognising that a 
quantum parameter estimation protocol can be written as a 
quantum network acting on a set of quantum systems, with 
repeated "black-box" couplings of the network to the system 
we wish to probe for the parameter (p [11]. A black-box is 
a function that can be univariate or multi-variate. When the 



function is multi-variate, a query to the black-box must take 
the form of multiple input parameters. Likewise, when the 
operator that describes the fundamental "atomic" interaction 
between the queried system and the probe is a two-body inter- 
action on the probe, then a query can consist only of precisely 
two input bodies. The scaling of the error in <p is then deter- 
mined by the query complexity of the network. The number 
of queries Q is not always identical to the number of physical 
systems N in the network. 

In Fig. lb-d) we consider three examples. The quantum 
network with univariate black-boxes in b) was analysed by 
Giovannetti, Lloyd, and Maccone [11]. Suppose that each 
grey box in Fig. 1 is a unitary gate Oj((p) = exp(—i(pHj), 
where j = 1,... ,N denotes the system, and Hj is a positive 
Hermitian operator. It is convenient to define the generator of 
the joint queries as J^glm = HjHj, because all Hj commute 
with each other. The number of queries Q is then equal to the 
number of terms in ,J#glm, or Q = N. In Fig. lc) the black- 
box is bi-variate. This is a type of Hamiltonian considered by 
Boixo, Flammia, Caves, and Geremia [9], and takes the form 
J^bfcg = LitLi Ey=i ®.Hfc. A physical query to a black- 
box characterized by Ojk = ex$(—i(pHj®Hk) must consist of 
two systems, labeled j and k. Since each pair interaction is a 
single query, the total number of queries is (^) = ^N(N — 1 ). 
Finally, in Fig. Id) we depict the network corresponding to the 
protocol by Roy and Braunstein [10]. It is easy to see that the 
number of terms in the corresponding generator J%rb is given 
by 2 N — 1, and the number of queries is therefore Q = 2 N — 1. 

A similar argument can be made to find the correct number 
of queries for all types of networks. The key principle is that a 
physical query to a quantum system consists of probe-systems 
that together undergo an operation, which can potentially en- 
tangle them. The entangling power of the black-box operation 
over multiple input systems accounts for the super-linear scal- 
ing of Q with N. Only when J%? does not have any entangling 
power across the input, we are guaranteed to have Q = O(N). 
This is in agreement with Refs. [9] and [10] where \/f (<p) 
scales super-linearly in N, but is always linear in Q, as de- 
fined here. Since we have a systematic method for increasing 
N (and Q) given the atomic interaction Hj, this uniquely de- 
fines an asymptotic query complexity of the network. Since 
both T and Q count the number of queries, this allows us to 
meaningfully compare the SQL with the Heisenberg limit. 

Given that in Eq. (2) ^F{(p) < Q, we have to find a gen- 
eral procedure that bounds Q, based on the physical descrip- 
tion of the estimation protocol in Fig. la). Previously, we 
showed that Q is the number of black-box terms in Jrif, and 
a straightforward choice for the resource count is therefore 
\(J^)\ < 0(Q). An important subtlety occurs when cor- 
responds to a proper Hamiltonian. The origin of the energy 
scale has no physical meaning, and the actual value of | ( Jtf? ) | 
can be changed arbitrarily. Hence, we must fix the scale such 
that the ground state has zero energy (equivalently, we may 
choose | {M 1 — hraaJ) \ , where h m m is the smallest eigenvalue, 
and / the identity operator). In most cases, this is an intuitive 
choice. For example, it is natural to associate zero energy to 
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the vacuum state, and add the corresponding amount of en- 
ergy for each added photon. Technically, this corresponds to 
the normal ordering of the Hamiltonian of the radiation field 
in order to remove the infinite vacuum energy. Slightly less 
intuitive is that the average energy of N spins in a GHZ state 
(I t) 8iV + I l) m ) /V2 is no longer taken to be zero, but rather 
N/2 times the energy splitting between | f) and | \). 

While the expectation value of Jtf? is easy to calculate, 
it is not the only way to obtain a bound of O(Q) from 
Other seemingly natural choices are the variance and 
the semi-norm. For example, if we write = E?Aj, the 

variance is (AJf) 2 = ((if A,) 2 ) - (LjAj) 2 = lf(Lj) - 
k (Aj) (Af,.) < cQ 2 for some positive number c and positive 
operator Lj. This gives AJ4? < O(Q), where e.g., in Ref. [9] 
Q = 0(N 2 ). Similarly, \{Jt?)\ = Zf\( A j)\ < 0(Q) since all 
expectation values are positive and finite. In other words, 
in terms of the scaling behaviour with Q, we can use either 
the variance or the expectation value. However, there are im- 
portant classes of quantum systems for which the variance of 
the energy diverges, such as systems with a Breit-Wigner (or 
Lorentzian) spectrum [12, 13]. Similarly, the semi-norm does 
not exist for a large class of states, such as optical Gaussian 
states. In these cases the resource count, and by implica- 
tion the scaling of the error, would be ill-defined. By con- 
trast, \(J^)\ always exists and is always positive. Also, from 
a physical perspective the higher-order moments do not de- 
scribe "amounts" in the same way as the first moment does, 
and refer instead to the shape of the distribution. This is a fur- 
ther argument that | (Jt?) | is the natural choice for the resource 
count. Sometimes, it is unclear how the query complexity is 
defined, for example when the estimation procedure does not 
involve repeated applications of the gates Oj((p), or when an 
indeterminate number of identical particles, such as photons, 
are involved. Nevertheless, the generator is always well- 
defined in any estimation procedure, and we can always use 
its expectation value to define the relevant resource count. 

The resource count in terms of | (Jff ) | is completely gen- 
eral for all possible quantum networks. For interactions U((p) 
where we include feed-forward and arbitrary unitary gates be- 
tween queries, we can use an argument by Giovannetti et al. 
[11] to show that = \{(dU(<p)/d<p)U^(<p))\ is unaf- 

fected by the intermediate unitary gates, and the scaling is 
therefore still determined by Q. 

After establishing the appropriate resource count, we are fi- 
nally in a position to prove the optimality of the Heisenberg 
limit for quantum parameter estimation in its most general 
form. The Fisher information can be related to a statistical 
distance s on the probability simplex spanned by p(x\<p). Con- 
sider two probability distributions p(x\(p) and p(x\(p) +dp(x). 
The infinitesimal statistical distance between these distribu- 
tions is given by ds 2 = f dx[dp(x\(p)] 2 / p(x\(p) [14, 15]. Di- 
viding both sides by (d(p) 2 we obtain 
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which relates the Fisher information to the rate of change of 
the statistical distance (i.e., the speed of dynamical evolution). 

When we count the resources that are used in a parame- 
ter estimation procedure, we must make sure that we do not 
leave anything out, and this can be guaranteed by including in 
our description the environment that the estimation procedure 
couples to. This reduces the quantum states to pure states, 
which means that we can use Wootters' distance [15] between 
quantum states as the statistical distance: 



*(V,0) =arccos(|(y|0>|), 



(4) 



where \\jf) and \<j>) are two pure states in the larger Hilbert 
space, and s(y,<j>) is the angle between them. The distance 
between the probe state p(0) and the evolved state p(<p) can 
then be represented by the pure states \y(0)) and re- 
spectively, and the unitary evolution is given by 



\\lf(( P ))=exp(-i( P M')\y(0)), 



(5) 



Here, we place no restriction on M 1 , other than fixing the en- 
ergy scale if necessary. We can place an upper bound on the 
derivative of Wootters' distance by evaluating the differential 
of s in Eq. (4) and using the Schrtidinger equation implicit in 
Eq. (5) [16]: 



£<K*>|. 



(6) 



Combining this with Eq. (3) and Eq. (2) leads to the Cramer- 
Rao bound 



1 / ds 
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T\dq>J -T|(jr)p 



(7) 



When all resources are used in a single-shot (T = 1) experi- 
ment, the error in <p is bounded by 



<S<p > 



1 



1(^)1 ' 



(8) 



Since | (Jf?) | is the resource count in the parameter estimation 
procedure, this is the Heisenberg limit. It is always positive 
and finite, and in the limit where | (Jf?) | — > there are no re- 
sources available to estimate <p, and 8q> cannot be bounded. In 
general, the bound is not tight. Indeed, only carefully chosen 
entangled systems can achieve this bound [11]. This com- 
pletes the proof of the optimality of the Heisenberg limit in 
the most general case. 

In addition to Eqs. (3) and (6), the Fisher information is also 
bounded by the variance of J4? according to F((p) < 4(AJf) 2 
[ 17] . This leads to a (single-shot) quantum Cramer-Rao bound 



<5<p > 



1 



(9) 



However, since AJf is not a resource count, such as the aver- 
age photon number, but rather a variance (or uncertainty) this 
is not the Heisenberg limit. In fact, it is Heisenberg's uncer- 
tainty relation for the parameter <p and its conjugate operator 
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J^ 7 . Any parameter estimation procedure must respect both 
bounds, and the Heisenberg limit in Eq. (8) may not be at- 
tained for a particular input state because the bound in Eq. (9) 
prevents it from doing so. 

The term "Heisenberg limit" was introduced by Holland 
and Burnett [7], who referred to the number-phase uncertainty 
relation in Heitler [18]. However, as our optimality proof and 
the subsequent discussion indicate, the Heisenberg limit is not 
an uncertainty relation, since it relates the uncertainty of the 
parameter to the first moment of the conjugate observable ,3^, 
rather than the second. It turns out instead that the Heisenberg 
limit is intimately connected to the Margolus-Levitin bound 
on the time it takes for a quantum system to evolve to an or- 
thogonal state [16, 19, 20]. To see this, we can formally solve 
Eq. (6) by separation of variables, yielding 

f<P , 1 r*l 2 % 1 

k d,f> -W)\l ds ri iWTi (10) 

We can therefore identify the Heisenberg limit with the 
Margolus-Levitin bound on the speed of dynamical evolution 
of quantum systems when Jtff is the Hamiltonian. The (gen- 
eralized) uncertainty relation, on the other hand, can be iden- 
tified with the Mandelstam-Tamm bound [16]. Both limits 
are completely general (assuming the existence of AJt?) and 
complement each other. 

Finally, we demonstrate that our proof applies to continu- 
ous variable systems as well as discrete systems, by consid- 
ering the procedure of Beltran and Luis [8]. The construc- 
tion is as follows: The evolution 0((p) is generated by an 
optical nonlinearity proportional to the square of the photon 
number operator n 2 acting on a single-mode coherent state 
1 1^(0)) = \a). The evolved state before detection is given by 
\\l/((p)) = exp(— /<pn 2 )|a), and the mean square error in (p is 
calculated as 8(p ~ \{n)~^l 2 = ||o:|~ 3 , to leading order in the 
average photon number («) . Since here the average energy is 
directly proportional to the average photon number, this pro- 
cedure seems to surpass the Heisenberg limit. To resolve this 
paradox, we note that the generator of translations in q> is not 
the photon number operator ft, but rather the higher-order non- 
linearity ,¥£ = ft 2 . The appropriate resource count is there- 
fore |(JT}| = (n 2 ), instead of the average photon number (n) . 
It is easily verified that to leading order 8(p is theoretically 
bounded by 1 / (n 2 ) = l/|a| 4 . Hence the parameter estimation 
procedure not only does not beat the Heisenberg limit, it does 
not attain it. 

Formally, we can attain the Heisenberg limit in this setup 
with the following modification of the input state and the 
measurement. Consider the single-mode input state \\j/o) = 
(|0) + \N)) /V2, where |0) denotes no photons, and \N) de- 
notes N photons. The state of the probe before detec- 
tion is then given by |v(<P)) = exp( — /<pn 2 ) 1 1//(0)) = (|0) + 
e -i<pN- |jv))/ sf2. We define the measurement observable X = 
\0)(N\ + \N)(0\. Hence, for the final state |y(<p)) we calculate 
(X) = (Y(p\X\y 9 ) =cos(N 2 (p) and AX = sin(N 2 (p). Using the 
standard expression for the mean squared error, we find that 

AX 1 

8<? = W)JM = ^- (U) 



Since \(Jf?)\ = (n 2 ) = jN 2 , this attains the Heisenberg limit. 
This is a formal demonstration that the Heisenberg limit can 
be attained according to quantum mechanics, even though we 
currently do not know how to implement it. 



In conclusion, we demonstrated that the Heisenberg limit 
is optimal for all parameter estimation procedures in quantum 
metrology, but it requires careful consideration as to which re- 
source is appropriate for expressing the scaling behaviour of 
the mean square error. The correct resource to take into ac- 
count is (the expectation value of) the generator of the trans- 
lations in the parameter. In the case of most optical phase 
estimation protocols this reduces to the average photon num- 
ber. Contrary to the origin of the term "Heisenberg limit", it 
is not a generalised uncertainty relation, but rather an expres- 
sion of the Margolus-Levitin bound on the speed of dynamical 
evolution for quantum states. 
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